ADP-05-08/T618, JLAB-THY-05-322 



Nucleon quark distributions in a covariant quark-diquark model 
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Spin-dependent and spin-independent quark light-cone momentum distributions and structure 
functions are calculated for the nucleon. We utilize a modified Nambu— Jona-Lasinio model in 
which confinement is simulated by eliminating unphysical thresholds for nucleon decay into quarks. 
The nucleon bound state is obtained by solving the Faddeev equation in the quark-diquark ap- 
proximation, where both scalar and axial-vector diquark channels are included. We find excellent 
agreement between our model results and empirical data. 
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I. INTRODUCTION 

The discovery in the late 1980's by the European Muon 
Collaboration (EMC) that the fraction of the spin of 
the proton carried by the quarks is unexpectedly small 
0, caused much excitement in the nuclear and particle 
physics communities. The "proton spin crisis" prompted 
many new experiments, leading to major new insights 
into the spin structure of the proton. However, a thor- 
ough theoretical understanding of the non-perturbativc 
parton distributions still remains an exciting challenge. 

In this paper we calculate the spin-independent 
and spin-dependent quark distributions in the frame- 
work developed by Bentz, Thomas and collaborators, 
in which proper-time regularization is applied to the 
Nambu— Jona-Lasinio (NJL) model Q in order to simu- 
late the effects of confinement . This model is attrac- 
tive because of its covariance and the transparent descrip- 
tion of spontaneous chiral symmetry breaking. We con- 
struct the nucleon as a bound state solution of the rela- 
tivistic Faddeev equation |j, la, IS L3 in the quark-diquark 
approximation , where both scalar and axial- vector di- 
quark channels are included. This quark-diquark descrip- 
tion of the single nucleon has the further advantage that 
it can be extended to finite baryon density [9j. We pay 
special attention to the polarized structure of the nu- 
cleon, comparing our results for the quark distributions 
with the empirical parameterizations. 



II. QUARK DISTRIBUTIONS 

The spin-dependent quark light-cone momentum dis- 
tribution in the nucleon, at leading twist, is defined by 
Eq. Q , where ip q is the quark field of flavour q, x is the 
Bjorken scaling variable and the subscript c reminds us 
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that only connected matrix elements are included. 



Afg(x) = P- I y^e 



( P , s \yj q (o) 1 + 15 Mr)\p,s) c . (i) 

We normalize the nucleon state vector according to non- 
covariant light-cone normalization: (p, s\^p q ^ + ip q \p, s) c — 
3. The spin-independent distribution, f q {x), is defined 
by the replacement 7+75 — > 7 + in Eq. (T|. To determine 
the quark distributions in this model, it is convenient to 
express Eq. in the form [HGIl 
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Tr 7 + 75 M g (p,fc) , (2) 



where M q (p, k) is the quark two-point function in the 
bound nucleon. Hence, within any model that describes 
the nucleon as a bound state of quarks, the distribution 
functions can be associated with a straightforward Feyn- 
man diagram calculation. 

The Feynman diagrams considered here are given in 
Fig. ^ where in our model the resulting distributions 
have no support for negative x. Therefore this is essen- 
tially a valence quark picture. By separating the isospin 
factors, the spin-dependent u and d distributions in the 
proton can be expressed as 

Au v (x) = Af° q/N (x) + \ Af° {D)/N (x) + 1 AQ N (x) 

+ l^ D)/ N(x) + ^=Af^ D)/N (x), (3) 



Ad v (x) = l - A.Q D)/N (x) + | Af a q/N {x) 

+ \Ar q{D)/N {x)-^j=AQ D]/N {x). (4) 



The superscripts s, a and m refer to the scalar, axial- 
vector or mixing terms, respectively, the subscript q/N 
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FIG. 1: Feynman diagrams representing the quark distribu- 
tions in the nucleon, needed in the evaluation of Eq. Q. The 
single line represents the quark propagator and the double line 
the diquark t-matrix. The shaded oval denotes the quark- 
diquark vertex function and the operator insertion has the 

form 7 + 7s 8^x — i (1 ± t z ) for the spin-dependent dis- 
tribution and 7 + 7s — + 7 + for the spin-independent case. 



implies a quark diagram and similarly q(D)/N a diquark 
diagram. Because the scalar diquark has spin zero, we 
have Af^ D y N (x) = and hence the polarization of the 
d-quark arises exclusively from the axial-vector and the 
mixing terms. Similar expressions hold for the spin- 
independent distributions, but in that case there is no 
mixing contribution (f% D y N = 0) Un- 
importantly, in this covariant framework, the Ward 
identities corresponding to number and momentum con- 
servation are satisfied, guaranteeing the validity of the 
baryon number and momentum sum rules 0, H3 ■ 



III. THE NUCLEON IN THE NJL MODEL 

The NJL model is a chiral effective quark theory that 
is characterized by a 4-Fermi contact interaction of the 

form, C[ — Gi ("0 Ti V>) , where the Fj represent ma- 
trices in Dirac, colour and flavour space and G{ are cou- 
pling constants 0- Applying Fierz transformations, the 
interaction Lagrangian can be decomposed into various 
interacting qq and qq channels. Writing only those terms 
relevant to this discussion, we have 



C = ijj (i$ — to) tp + C Ii7r + Ci. s + Ci. a 



(5) 



where m is the current quark mass. The interaction 
terms are given by 

=\G, ((?V) 2 - (^7sr ^) 2 ) , (6) 
Ci, s = G s (^ l5 Cr 2 f3 A i?) (V C- X l5 r-20 A i>) , (7) 
Hi, a = G a (^^CnT 2 p A ^ T ) (Vc-V^t^t/;), (8) 



where /3 — y | A A (A = 2, 5, 7) arc the colour 3 matrices 
and C = ?727o- The familiar term Ci^ generates the 
constituent quark mass, M, via the gap equation and 
the pion as a qq bound state. The terms £/ iS and Ci^ a 
represent the interactions in the scalar (J 11 — Q + ,T — 
0, colour 3) and axial-vector (J" = 1 + ,T = 1, colour 3) 
diquark channels and are used to construct the nucleon as 



a quark-diquark bound state. The couplings G^, G s and 
G a arc related to the original couplings, Gi, via the Fierz 
transformation, but we use them here as free parameters 
which will be fixed by the properties of the pion and the 
nucleon. 

Solving the appropriate Bethe-Salpeter equations, the 
standard NJL results for the diquark i-matrices are ob- 
tained 0, As explained in Ref. @, these can be 
accurately approximated by the forms 



T s (q) = AiG s 



Ws 



ql-Mf 



r^(q)=4dG a g^ 



A/ 2 



(9) 
(10) 



which we also use here. The masses of the diquarks 
M s , M a and their couplings to the quarks g s , g a are de- 
fined as the poles and residues of the appropriate full 
diquark i-matrices. 

The nucleon (quark-diquark) t-matrix satisfies the 
Faddeev equation 



T=Z + Zn N T = Z + TI1 N Z, 



(11) 



where Z is the quark exchange kernel and Hn the product 
of a quark propagator and a diquark t-matrix. In the 
non-covariant light-cone normalization used already in 
Eq. the quark-diquark vertex function, Tat, is defined 
by the behaviour of T near the pole 



T 



p + ^ep r^v r^v 



p+ 



(12) 



where e p = is the light-cone energy. Substituting 
this result into Eq. Ijllfl gives the homogeneous Faddeev 
equations for the vertex functions 



N 



zn N v 



N 1 JVj 



and T N = T N U N Z 



(13) 



For this investigation we restrict ourselves to the static 
approximation, where we neglect the momentum depen- 
dence of the quark exchange kernel, Z. Including both 
scalar and axial- vector diquark channels, Z takes the fol- 
lowing form in the colour singlet and isospin-i channel: 



z = _3_ ( 1 v / 37/.'75 
M V\/3~757 M "7m'7 m 



(14) 



The quantity LTjv effectively becomes the quark-diquark 
bubble graph: 



n N (p) = J 



where 



r{q) 



Ts{q) 



-(p-k)S(k), 







tT(?) 



(15) 



(16) 
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The eigenfunction of the kernel K = Z Hn, in Eq. (|13|l . 
has the following form, up to normalization: 



r(p, s) 



a i 



oil 



Mr 



■75 +a3 7 AI 75 



un(p,s), (17) 



where the upper and lower component refer to the scalar 
and axial- vector diquark channels, respectively and un 
is a free Dirac spinor with mass Mm- We choose the nor- 
malization umUm = 1 = IT. 1 Inserting this form into 
Eq. (fT3*|l gives three homogeneous equations for the a's 
and the nucleon mass Mjv is determined by the require- 
ment that the eigenvalue of if, in Eq. I|13[) . equal 1. 

The normalization of the vertex function follows from 
the definition given in Eq. 1)12(1 . we obtain 



Tn(p, s) 



p- 



where 



P- 



m n r( p )^Mr(p) 



(18) 



(19) 



As with any non-renormalizable theory a regulariza- 
tion prescription must be specified to fully define the 
model. We choose the proper-time regularization scheme 
IE El El , where loop integrals of products of prop- 
agators are evaluated by introducing Feynman parame- 
ters, Wick rotating and making the denominator replace- 
ment 



1 



1 



{n-l)\ 



l/(Am) 2 
l/(Ac/v) 2 



dTT n - x e- rX , 



(20) 



where Aju and A;yy are, respectively, ultraviolet and in- 
frared cutoffs. The former has the effect of eliminat- 
ing unphysical thresholds for hadron decay into quarks, 
hence simulating an important aspect of confinement 

M. 



IV. RESULTS 

The parameters of the model are Am, A^y, m, G„, 
G s and G a . The infrared scale is expected to be of or- 
der Aqcd and we set it to Ajr = 0.28 GcV. This is 
slightly larger than our previous work [3j, because our 
studies of the saturation properties of nuclear matter 
favour this 0. The parameters m, Auv and are 
determined by requiring M = 400 MeV via the gap equa- 
tion, f„ = 93 MeV from the familiar one loop pion decay 
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FIG. 2: Spin-independent valence u and d distributions mul- 
tiplied by Bjorken x. There are three curves for each quark 
flavour, with the lower curve of each type representing the d 
distribution. The dotted line is the model prediction at the 
NJL scale of Q 2 , = 0.16 GeV 2 and the solid line is the result 
after QCD evolution to the scale Q 2 = 5.0 GeV 2 . The dashed 
line is the empirical parametrization of Ref . [Tsll , at the scale 
Q 2 = 5.0 GeV 2 . 



diagram and m n — 140 MeV from the pole of the qq t- 
matrix in the pion channel. This gives m — 15.3 MeV, 
A uv = 0.66 GeV and G^ = 17.81 GeV~ 2 . The cou- 
plings G s and G a are determined by reproducing the nu- 
cleon mass M N = 940 MeV as the solution of Eq. 1(13)1 
and satisfying the Bjorken sum rule within our model, 
where g A = 1-267. We obtain G s = 8.41 GeV" 2 and 
G a = 1-36 GeV -2 . With these model parameters the di- 
quark masses are M s = 0.65 GeV and M a = 1.2 GeV and 
the coefficients in the nucleon vertex function, Ea. (|17ll . 
are (a l7 a 2 ,a 3 ) = (-0.35,-0.0088,0.47). 

To compare the predictions of the model with experi- 
mental data as well as the empirical parameterizations, 
it is necessary to determine the model scale, Qg. We 
do this by optimizing Qq such that the spin-independent 
distribution, u v (x), best reproduces the empirical param- 
eterization after Q 2 evolution. We find a model scale of 
Qq = 0.16 GeV 2 which is typical of valence dominated 
models 0E[l7i. 

Results for the spin-independent and spin-dependent 
valence u and d distributions are presented in Figs. 
and |3 respectively. We show the predictions at the 
model scale and after QCD evolution 2 to Q 2 = 5 GeV 2 , 
where they are compared to empirical parameterizations. 
We find excellent agreement between the model results 
and the parameterizations. In comparison with the pure 



The conjugate vertex function, T, which is a left eigenfunction 
of K = IIjv Z, is obtained by taking the ordinary hermitian 
conjugate of V and introducing a minus sign for the axial- vector 
components. 



2 We utilize the computer program of Ref. |2(1 for the spin- 
independent case and of Ref. |2l| for the spin-dependent case. 
We choose DGLAP evolution with N f = 3, Aqcd = 250 MeV 
in the MS renormalization scheme up to NLO. 
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FIG. 3: Spin-dependent valence u and d distributions multi- 
plied by Bjorken x. The curves are as in Fig- [21 with empirical 
parameterizations taken from Ref. |l9(|. 



scalar model 0, 12^, the agreement has improved sub- 
stantially, especially for the spin-dependent case. 

Our model results for the first polarized moments are 
Au v = 0.924 and A d v = —0.343 which agree with the 
values Au v = 0.926 ± 0.014 and Ad v = -0.341 ± 0.018 
determined from the axial coupling constants of octet 
baryons discussed in Ref. j23|. This emphasizes the im- 
portance of including axial- vector diquark correlations, 
since the pure scalar model would give a vanishing A d v 
and a somewhat smaller A u v . The spin sum in our model 
is A§ = 0.581, which is smaller than the result of the 
pure scalar model, but still somewhat larger than the ac- 
cepted value of A§ = 0.213±0.138 This discrepancy 
primarily reflects the absence of the ?7(1) axial anomaly 

The behaviour of structure function and hence quark 
distribution ratios at large x has been an area of consid- 
erable debate |26L and is one of the regions where per- 
turbative QCD (pQCD) offers firm predictions |2g. Ex- 
perimentally, the ratio d(x)/u(x) is surprisingly poorly 
known j2^. In the limit x — > 1 it is thought to lie 
somewhere between 0, the prediction based on scalar 
diquark dominance [3(j and |, the pQCD result |28j . 
Analysis in Ref. [2(j favours the pQCD prediction. The 
same predictions also hold for the spin-dependent ratio, 
Ad(x)/Au(x), as x approaches 1. 

In Fig. 0] we plot our results for the ratios d v (x)/u v (x) 
and Ad v (x) / Au v (x), together with the ratios of the em- 
pirical distributions. The x — > 1 limit of the spin- 
independent ratio is in agreement with the pQCD result. 
The spin-dependent ratio, however, approaches ~ — j^, 
the opposite sign to the pQCD prediction. Although the 
empirical parameterizations are constrained to give for 
these ratios as x — > 1, we note that the systematic er- 
rors in both empirical ratios are very large in the region 

x > 0.5 [UlilHIal- 

It is important to note that the pQCD predictions 
for the mixed flavour ratios are somewhat model depen- 



FIG. 4: Mixed flavour ratios for spin-independent and spin- 
dependent distributions. There are two curves for each ratio, 
with the lower curves the polarized result. The curves are 
as in Fig. [21 with spin-independent parameterizations taken 
from Ref. |1S| and the spin-dependent from Ref. [is| 
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FIG. 5: Single flavour ratios (Aq + Aq) / (q + q) where q G 
(if, d), at the scale Q 2 = 5.0 GeV 2 . The lower curve represents 
the d-quark distribution. The experimental results are from 
Hall A at Jefferson Lab |27| (solid squares) and Hermes pjUl 
(solid stars). 



dent, as assumptions have to be made about the relative 
strengths of the u and d contributions to the nucleon 
wavefunction. A more rigorous pQCD prediction, rely- 
ing only on helicity conservation, is possible for the single 
flavour ratios Au(x)/u(x) and Ad(x)/d(x). Perturbative 
QCD predicts that both these ratios should approach 1 
for large x, which would require a change of sign in the 
Ad distribution. 

In Fig. [S] we plot our results for the ratios 
(Aq + Aq) / (q + q) where q € (it, d). Since we wish 
to compare these ratios directly to recent experimen- 
tal data, we include sea quark distributions generated 
through the Q 2 evolution. In the x — > 1 limit our model 
ratios approach w 0.8 for the M-quark and w —0.25 for the 
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FIG. 6: Polarized structure functions gi p and g\ n at Q — 
5GeV 2 . The solid line is the model prediction, with the 
lower curve corresponding to g\ n . The shaded areas rep- 
resent the empirical parameterizations with uncertainties of 
Ref. [3^1 . at the same scale. The experimental data, with 
1 < Q 2 < lOGeV 2 , is from SMC Jp (open stars), SLAC 
E143 |23 (open circles) and JLab |27fl (solid squares). 



FIG. 7: Structure function ratios A\ p and Ai n , at the scale 
Q 2 — 5 ,GeV 2 . The solid line is the model prediction, with 
the lower line corresponding to Ai„. The experimental points 
represent data with 1 < Q 2 < 10 GeV 2 from Hermes |37)| 
(closed star) SMC [3^| (open star), E143 (open circles) 
and JLab 27] (solid squares). We do not plot the E155 
results as they are similar to those from E143. 



d-quark. This seeming contradiction to pQCD has also 
been suggested by recent experiments by the Jefferson 
Lab Hall A collaboration |27J, |35j , with our predictions 
consistent with their experimental results. This data is 
also shown in Fig. [SJ 

In Fig. we give our results for the spin-dependent 
structure functions g\ p (x) and g\ n {x). The parameteri- 
zations of Ref. are also included as the shaded ar- 
eas, which indicate the empirical uncertainties. Our re- 
sults compare well with the empirical parameterizations, 
agreeing within uncertainties for the region x > 0.25. 
Comparison with experiment is also favorable, although 
the experimental determination for gi n {x) is less certain. 

Model results for the asymmetries Ai p (x) and Ai n {x) 
are given in Fig. \7\ where we see good agreement in the 
valence region. However, for small x, A\ p (x) is slightly 
too large, because of the enhancement of gi p (x) in the 
same region. This is most likely associated with the omis- 
sion of the effects of the axial anomaly in the present 
work. It is also clear from the experimental data that 
the uncertainties in these ratios at large x, are still sig- 
nificant. 

V. CONCLUSION 

Using a covariant quark-diquark model for the nucleon, 
including both scalar and axial- vector diquark channels, 
we calculated the spin-independent and spin-dependent 
quark light-cone momentum distributions and structure 
functions. A key feature of the framework is that it pro- 
duces quark distributions that have the correct support 
and obey the number and momentum sum rules. The 
model also incorporates important aspects of confine- 
ment by eliminating unphysical thresholds for nucleon 



decay into quarks. 

Highlights of our results are obtaining values for the 
polarized first moments of the quark distributions Au v = 
0.924 and Ad v — —0.343, in agreement with those ob- 
tained from axial couplings of octet baryons. We also 
obtain excellent agreement with empirical parameteriza- 
tions of the valence quark distributions. We paid special 
attention to the single flavour ratios ( Aq + Aq) / (q + q) 
and the asymmetries A\ p and A\ n , finding good agree- 
ment with recent experimental results from JLab. 

These results indicate that diquark correlations are an 
essential feature of the non-perturbative structure of the 
nucleon. In particular, the admixture of axial-vector di- 
quarks, though small, is essential to obtain the observed 
agreement with empirical data. 

Finally, we would like to mention that a very important 
advantage of this covariant quark-diquark model is that 
it can be readily extended to the case of finite nucleon 
density. The results presented in this paper strongly sug- 
gest that this model should provide a reliable basis from 
which to begin investigation of the medium modifications 
of both spin-independent and spin-dependent structure 
functions. 
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